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Abstract 

The phase diagram of five-dimensional SU(2) gauge theories is explored 
using Monte Carlo simulations of the theory discretized on a Euclidean lat- 
tice using the Wilson plaquette action and periodic boundary conditions. 
We simulate anisotropic gauge couplings which correspond to different lat- 
tice spacings 04 in the four dimensions and 05 along the extra dimension. In 
particular we study the case where 05 > 04. We identify a line of first order 
phase transitions which separate the confined from the deconfined phase. 
We perform simulations in large volume at the bulk phase transition staying 
in the confined vacuum. The static potential measured in the hyperplanes 
orthogonal to the extra dimension hint at dimensional reduction. We also 
locate and analyze second order phase transitions related to breaking of the 
center along one direction. 
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1 Introduction 



Our interest in studying five-dimensional gauge theories comes from extensions 
of the Standard Model called Gauge-Higgs unification. The idea is to identify 
the Higgs boson with (some of) the extra dimensional components of the gauge 
field. Since five- dimensional gauge theories are non-renormalizabl^ an ultra-violet 
cut-off is mandatory. The lattice regularization provides such a gauge-invariant 
cut-off (the inverse lattice spacing) and is therefore a natural setup to study these 
theories. 

Here we consider the five-dimensional pure SU(2) gauge theory and investi- 
gate its non-perturbative phase diagram through Monte Carlo simulations. The 
theory is discretized on a Euclidean space-time lattice. The seminal work done by 
Creutz [ilj using the Wilson plaquette gauge action [2] demonstrated the existence 
of (at least) two bulk phases in the non-perturbative phase diagram of the theory 
in infinite volume. There is a confined phase at large values of the gauge coupling 
and a deconfined phase when the gauge coupling is small. 

The question we are after is where dimensional reduction from five to four 
dimensions occurs. Possible mechanisms, which have been proposed are compacti- 
fication or localization. In compactification models the size of the extra dimension 
is small and the system is described by an effective four- dimensional low energy 
theory valid for energies much below the inverse compactification radius. The ex- 
istence of a minimal physical size of the extra dimension, below which the theory 
becomes four- dimensional has been suggested by lattice simulations in |3]. The 
extra dimension can even "disappear" when the correlation length of the four- 
dimensional theory grows exponentially fast with the size of the extra dimension. 
This is the mechanism of D-theory which relies on the existence of the 

Coulomb phase and has been recently studied on the lattice J6j. 



The most prominent of localization models in fiat spacci is the domain wall 
construction for fermions which has been shown to be equivalent to an orb- 
ifold construction [8]. A possible localization mechanism for gauge field has been 
proposed in |[9] and relies on a mechanism that confines the theory in the bulk 
and deconfines it on a boundary (or domain wall) of the extra dimension. This 
mechanism has been studied on the lattice in [10], where it was found that indeed 
localization occurs but the low energy effective theory contains not just the local- 
ized zero-modes but also higher Kaluza-Klein modes. Similarly localization can 
can also be realized if a layered phase exists [TT|IT^. 

The phenomenological signature of extra dimensional models depends on the 

^ Equivalently one can say that five-dimensional gauge theories are triviaL 
^ We do not discuss here the case of warping, which modifies the metric. 
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mechanism of dimensional reduction. In the compactification case a striking evi- 
dence would be the appearance of Kaluza-Klein modes. For example, so far the 
non-observation of an excitation Z' of the Z boson at LHC [T3|[H] has put a lower 
bound on their mass of about 1 TeV (assuming it has the same couplings as the 
Standard Model particle). This translates into a bound on the inverse compact- 
ification radius. In the localization case, there is a spectrum of four- dimensional 
localized modes at low energies with a mass gap, which is set by the domain 
wall height M d, separating the higher modes. The latter can be localized modes 
forming like a Kaluza-Klein tower or bulk modes [8|[T0]. 

Recently [T5l[T6] the phase diagram of the five- dimensional pure SU(2) gauge 
theory on the torus discretized with the Wilson plaquette action and anisotropic 
gauge couplings has been investigated in the mean-field approximation including 
corrections [T71[T8]. The mean-field phase diagram exhibits three phases, besides 
the confined and deconfined phases there is a layered phase when when the lattice 
spacing along the extra dimension 05 is larger than the one in the usual four 
dimensions 04. It is possible to take the continuum limit approaching a critical 
line (a line of second order phase transition^ ) at the phase boundary between 
the deconfined and the layered phase. The lattice spacings go to zero by keeping 
the anisotropy and the ratio of the vector to the scalar masses fixed. In this limit 
the hyperplanes orthogonal to the extra dimension decouple from each other and 
the theory turns four- dimensional. In the present paper we show what happens 
when the full theory is simulated by Monte Carlo methods. 

The paper starts with Section[2]containing definitions of the lattice theory and 
observables used in this study. In addition to observables (plaquettes and Polyakov 
loops and their susceptibilities. Binder cumulants) signaling phase transitions, we 
compute the static potential in the four- dimensional hyperplanes orthogonal to the 
extra dimension. Using the latter we define renormalized couplings from the static 
force and its derivative, which give us information about dimensional reduction. In 
Section |3] we discuss the case of isotropic couplings (where the lattice spacing is the 
same in all directions) as a starting point for the investigations of the anisotropic 
case in Section HI We study the part of the phase diagram where 05 > 04. We find 
first order bulk phase transitions (Section 14. ip . which separate the confined from 
the deconfined phase, and second order phase transitions related to breaking of the 
center symmetry along one lattice dimension (Section l4.2p . In Section we show 
our results for the static potential measured at two bulk phase transition points 

■^On the lattice M is given by the inverse lattice spacing. 

^ The existence of a second order phase transition for five-dimensional Yang-Mills theories 
has been hinted by the epsilon expansion [lUEO] but has been sofar elusive in lattice Monte 
Carlo simulations. 
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in the confined vacuum, which hint at dimensional reduction. Our conclusions 
and outlook are presented in Section [5l Appendix |A] summarizes the formulae for 
the 3-loop perturbative running of the coupling derived from the static force in 
the four- dimensional Yang-Mills theory, which we use to compare with our lattice 
data. 



2 Definitions 

We consider a Euclidean lattice in five dimensions. The dimensions are labelled 
by M = 0, 1, 2, 3, 5. The lattice size is V = x x L5, where Lt, Ls and L5 are 
the number of lattice points in the temporal (M = 0), spatial (M = 1,2,3) and 
extra (M = 5) dimensions respectively. We use Roman capital letters M, N, . . . 
to label all five dimensions and Greek letters yU, z/, . . . to label the four dimensions 
0, 1, 2, 3. The lattice points are denoted by a; = (xq, x, X5) and the gauge links by 
U{x, M). The Wilson lattice action for a SU(2) gauge theory is 



- J2 - U{x- /i, z/)} + 7 5^ tr{l - U{x- /i, 5)} 



7 



, (2.1^ 



where U{x] M, N) is the oriented plaquette at point x in directions M and 
and we use the property that the trace of SU(2) matrices is real. The parameter 
7 is the anisotropy parameter. Instead of 7) we will also use the equivalent 
parameter pair 

/54 = - , /95 = /57- (2.2) 

7 

The isotropic lattice corresponds to 7 = 1. If 7 7^ 1 there are two different lattice 
spacings, 04 in the /^-dimensions and 05 in the extra dimension. At the classical 
level the relation 7 = 04/05 holds. The relation between the lattice coupling (3 
and the bare dimensionful gauge coupling is [3] 

P = ^,(iV = 2). (2.3) 

The coupling is an effective coupling at the cut-off scale and an attempt to 
define a continuum limit 04 — )■ by keeping g^ constant (up to a slowly varying 
renormalization factor) leads to decreasing /3. Eventually the value Pc at the 
phase transition between the deconfined and confined phase is reached, meaning 
that the lattice spacing cannot be further reduced. The continuum limit can 
therefore only be taken at g^ = 0, i.e. the theory is trivial [21]. Triviality can 
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also be understood by considering the 1-loop renormalization of the effective four- 
dimensional coupling gl = gl/^L^a^) [221123]. 

In the Monte Carlo simulations of Eq. (12. ip we use the heatbath algorithm 
proposed by Fabricius and Haan [2l] and independently by Kennedy and Pendle- 
ton [25] in combination with the overrelaxation algorithm, cf. [26]. One updat- 
ing unit (or iteration) is typically defined as the sequence of Ls/2 overrelaxation 
sweeps followed by one heatbath sweep. In order to check the ergodicity of our 
simulations we have repeated some of them using Creutz's version of the heatbath 
algorithm [27] and a mixture of the Kennedy- Pendleton and Creutz one. The sim- 
ulation results obtained with these three variants of the algorithm agree very well. 
The statistical analysis of our simulations is usually done with the method of [2S] 
but in some cases we use the bootstrap method. 

The observables that we measure to investigate the phase diagram are the 
traces of the plaquettes, separately for the plaquettes orthogonal to and along the 
extra dimension 



Plaq4 = ^^^tr{f/(x;/i,z/)} , Plaqg = ^ ^ ^tr{f/(x; /i, 5)} 



and their susceptibilities 



XPiaq, = V ((Plaq, - (Plaq,))') , z = 4, 5 . 



(2.4) 



(2.5) 



The plaquette observables provide signals for bulk phase transitions. In order to 
study phase transitions related to the breaking of the center symmetry we also 
measure Polyakov loops. For example, the Polyakov loop along the time dimension 
and its susceptibility are defined by 



PolyT = Y' 



{LT-l)a4 

J2tT H U{x,0) 

X,Xc, X(} = 



V 2\ 

XPoiy^ = — ((Poly^ - (Poly^)) ) . 



(2.6) 



(2.7) 



Polyakov loops winding in other directions are similarly defined. We will also 
consider the temporal Polyakov loop without averaging over the extra dimension, 
explicitely 



{LT-l)a4 

^tr n U{x,0)\ 

X xo=0 



(2.8) 
(2.9) 
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It will be clear in the following why we use two different definitions Eq. (12 .Gp and 
Eq. (12.81) of the Polyakov loop. Other useful quantities to investigate the order 
of phase transitions are the Binder cumulants [29]. For example, the fourth order 
Binder cumulant of the temporal Polyakov loop is 

(Polyt) 

S4 = 1-^ (2.10) 

3 (Poly^)^ 

and similarly for other quantities. 

Wilson loops are traces of product of links along rectangular paths (we con- 
sider only on-axis Wilson loops). We measure Wilson loops W{r,T) of size T in 
the temporal direction and r along one of the spacial dimensions. On anisotropic 
lattices there are two types of loops, depending whether the distance r is taken 
along the extra dimension or the other spatial dimensions. In the measurement of 
the Wilson loops, the temporal links are replaced by their one-link integrals [30] 
and the spatial links are replaced by the links obtained after two levels of HYP 
smearing ^T] . The HYP smearing of the spatial links is restricted to staples ex- 
tending in the spatial directions, e.g. for a link in direction /x = 1 the staples are 
in directions M = 2, 3, 5. We use the HYP parameters 

ai=0.8, a2 = 0.56, 03 = 0.24 (level 1) , (2.11) 
ai = 0.92, ^2 = 0.68, 03 = 0.28 (level 2) , (2.12) 

optimized by maximizing the minimum plaquette. We find that these values do 
not strongly depend on the gauge couplings. We increase the levels of smearing up 
to four by using for level 3 and 4 the same parameters as for level 2. We checked 
numerically that the average and minimum plaquette increase with the number of 
smearing levels since our parameters ai are larger than the perturbative bound on 
the APE smearing parameter a < 0.75 [32j. From the Wilson loops we compute 
the effective masses 

a,m,s{r,t + a,/2) = J^^^'^^\ ^ '^^a,V{r) (2.13) 

\{W{r,t + ai)); 

and extract the static potential V{r) from a plateau average of the effective masses 
at large enough t, cf. [33]. In order to determine where to start the plateau average 
for each value of r we perform a fit 

meff(r,t + a4/2) = E + he'^^^^'^^ (2.14) 

with fit parameters b and A. The plateau average starts at the earliest time 
to when the exponential correction in Eq. fl2.14p is smaller than 1/4 of the the 
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statistical error of the effective mass. The plateau ends before the time t when 
either the difference of the effective mass at time t to the one at to is larger than 
the error of the effective mass at to ot the error of the effective mass at time t is 
larger than twice the error of the one at to. 

In our analysis we focus on the static potential along spatial directions or- 
thogonal to the extra dimension and we call this potential V{r). We define the 
static force through the symmetric derivative of the potential 

F(r-a4/2) = {V{r) - V{r - a4)}/a4 . (2.15) 

The scale ro/a4 [31] is determined through the numerical solution of 

r'F{r)l^^^ = 1.65, (2.16) 

which we find using a 2-point interpolation of the force F{r) = fo + f2/i^^- A 
comparison with a 3-point interpolation adding a f^/r^ term allows to estimate 
the systematic error. From the static force, a running coupling can be defined (in 
the so called qq-scheme) 

«qq(l/0 = ^r'F{r), (2.17) 

where Cp = 3/4 for gauge group SU(2). The perturbative expansion up to three 
loops of Oqq in four-dimensional Yang-Mills theory is summarized in Appendix 
lAl By taking a further derivative we compute the slope [35] 

c(r) = ir3F'(r), (2.18) 

which also defines a running coupling (at small distances). On the lattice we 
compute c(r) through 

c(r) = ^r^V{r + ai) + Vir-ai)-2Vir)}/al. (2.19) 

At large distances, the slope can be compared to the result of effective bosonic 
string theory [36l[37], which yields the asymptotic value 

= (2-2°) 

where D is the number of space-time dimensions. Thus c can serve to determine 
the number of dimensions. We remark that due to the factor in its definition, 
the noise-to-signal ratio of c(r) rapidly deteriorates when r grows. 

Since most of our large volume simulations are done at 7 < 1 in the confined 
phase, it turns out that the static potential along the extra dimension V^lr) is 
difficult to extract. The reason is that the values of a4V5(r) we measure at r = 205 
are typically 1 or larger and the dimensionless string tension (the coefficient of the 
term in V^lr) which is linear in r) a^a^a^ contains the large lattice spacing 05. 
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Figure 1: Plaquette and Polyakov expectation value as function of j3 on isotropic lat- 
tices. 



3 Isotropic couplings 

The aim of this section is to review known results for and highhght features of the 
isotropic case (7 = 1) that will be of help in the understanding of the anisotropic 
(7 7^ 1) case. We will start with the description of the system at zero temperature 
in the "infinite" volume limit. Afterwards we will analyze how any kind of finite 
size affects the results. 

The phase diagram of the five- dimensional isotropic SU(2) lattice gauge the- 
ory is well known since the early work of Creutz [1]. In the zero temperature 
setting and at infinite volume (both spatial and in the extra-dimension) the phase 
diagram of the model is characterized by the presence of a bulk first order phase 
transition at a critical value of the lattice coupling /3 = /3c ~ 1.64. The bulk phase 
transition is characterized by the appearance of a strong hysteresis region in the 
expectation value of the plaquette, see Fig. [H that survives in the infinite volume 
limit, signal of the presence of latent heat at the transition point. The two regions 
separated by /3c have a different behavior of the Wilson loop, with a string tension 
(the coefficient a of the linear term or in the static potential V{r)) different from 
zero in the confined phase for P < I3c and a string tension equal to zero in the 
deconfined phase above /3c. On a finite (but large) lattice this coincides with a 
distribution of the Polyakov loop (for all directions and not taking the absolute 
value) being always peaked in zero in the first case, while always showing a double 
peak structure in the second case. 

In Fig. [T]we report an example of the typical plots showing the behavior at the 
phase transition both for the expectation value of the plaquette and of the absolute 
value of the temporal Polyakov loop Poly^^ defined in Eq. (12. 6p . "Hot Start" or 
"Cold Start" means that the simulation starts from an initial configuration with 
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Figure 2: Expectation values of the plaquette Plaq4 and its susceptibility XPiaq4 as 
function of (3 on lattices with L5 = 2 and different values of Lj" = Lg = 6, . . . ,16. 

random gauge links or gauge links set to the identity, respectively. We perform 
10^ measurements separated by 10 updating iterations. 

3.1 One "small" dimension 

Now we investigate the effects of the finite lattice size on our observables and how 
they affect the bulk phase transition. As it is well known there is no sign of first 
order bulk phase transition in four-dimensional SU(2) Yang-Mills theory, instead 
a sharp crossover separates the weak and the strong coupling regime. It is then 
legitimate to ask which is the fate of the five-dimensional bulk phase transition 
when we consider the five-dimensional system compactified along one "small" 
direction. 

If we consider a geometry that has one small direction while all the other 
directions are sent to infinity, namely L5 -C Lt = Lg, we notice that the bulk phase 
transition at /3 = 1.64 "evaporates" and lets a smooth cross-over appear, cf. the 
early work of [SS]- This is shown by the behavior of the four-dimensional plaquette 
Plaq4 Eq. (12. 4p and of its susceptibility XPiaq4 Eq. (12. 5p in Fig. |2l for L5 = 2 and 
different choices of the four- dimensional lattice size Lt = Ls = 6, . . . , 16. There 
is no dependence nor on the initial condition (hot or cold start) of the simulation 
neither on the lattice volume. The cross-over region is close to (3 = 1.48 and will 
be investigated below. The susceptibility of the plaquette shows a peak at (3 values 
reminiscent of the bulk phase transition but whose magnitude is independent on 
the volume. We extend the simulations for our largest 2 x 16^ lattice up to /3 = 2.4 
but we do not see any sign of the cross-over of the four- dimensional SU(2) theory, 
which is at /3 ^ 2.3 
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Volume 


13, (Poly^) 


/3e (Polyr(l)) 


Xc (PoIyt) 


Xc (Polyr(l)) 


2 X 44 


1.50061(19) 


1.50369(23) 


11.86(4) 


3.842(15) 


2x6^ 


1.49097(13) 


1.49249(14) 


30.00(13) 


6.571(19) 


2x8^ 


1.48748(11) 


1.48839(10) 


56.9(4) 


9.44(5) 


2 X 10^ 


1.48570(10) 


1.48621(12) 


93.1(8) 


12.19(9) 


2 X 16^ 


1.483952(16) 


1.484170(18) 


257.8(1.6) 


20.88(10) 


2 X 00^ 


1.48280(1) 


1.48280(3) 







Table 1: For the center breaking phase transition at 7 = 1, Lx = 2, we Hst the critical 
values of the coupling /3c and the susceptibility of the temporal Polyakov loop Xc as 
function of the volume. We use two definitions of the Polyakov line, Eq. ()2.6p and 
Eq. ra. 

Our simulations hint at the existence of a minimal lattice siz^ 

2 < L'"''^(7 = 1) < 6 , (3.1) 

for which if one or more of the lattice dimensions are smaller then L™'"(7 = 1) we 
are unable to detect any sign of the bulk phase transition. 

Apart from the study of the plaquette, other interesting signals come from the 
study of the expectation value of the Polyakov loops and their susceptibility in the 
geometry with one small direction. It is interesting to notice that together with the 
disappearance of the bulk phase transition also the double peak structure of the 
Polyakov line distribution along the large directions disappears. On the other side 
we expect, and find, a proper center breaking structure along the small direction, 
which we choose for this study to be the temporal one. In order to avoid the 
appearance of the bulk phase transition we always request that < L™™(7 = 1). 
Notwithstanding the fact that we need to fix to a finite value, we can send the 
orthogonal space size L = Lg = L5 to infinity and obtain a system that undergoes a 
proper phase transition. We performed a study of the finite size scaling properties 
of this breaking by investigating both the distribution of the temporal Polyakov 
loop and its Binder cumulant (as discussed in [6]) as function of the orthogonal 
lattice size, while keeping = 2 fixed. 

We did a scan in /3 for the different choices of the lattice volume and we 
used a reweighting technique to obtain a more dense scan. The technique is the 
multi-histogram reweighting method [IQ]. Its main idea is to find iteratively the 
free energy of the system as function of the couplings using all the simulations 

^ We are exposing the dependence on 7 in order to extend later the definition of the minimal 
size to the case 7 7^ 1. 
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Figure 3: The susceptibility of Poly-p for the volume 2 x 16^. Comparison of the directly 
simulated data set (plusses) and the reweighted data (crosses) . 

performed in the region of interest. Once the procedure has converged, we have 
access via the free energy to all the physical quantities that can be normally 
extracted in a Monte Carlo simulation, for each value of the coupling constants 
in the region. Our implementation is able to handle simultaneous reweighting in 
both the couplings /34 and /Ss, even though in this section we are interested in the 
(3^ = 13^ case. 

In order to determine the error of our estimates we use the bootstrap method 
and do a reweighting analysis over each bootstrap sample. The underlying idea 
is to consider the estimates coming from each bootstrap sample as independent 
measurements and to use the spread of these measurements to define the error of 
the average of the estimates. In Fig. [3]we clearly see that our reweighted estimates 
for the susceptibility XPoiyj, (crosses, they appear like a band in the figure) Eq. (12.71) 
are more precise than the measurements themselves (plusses). This is however not 
surprising, since the reweighting technique uses for each estimate the statistical 
information coming from all the simulated points. 

The critical coupling I3c is defined as the coupling at which the susceptibility 
XPoiyy has its maximum Xc = XPoiyr(/3c)- The results for /3c{L) and Xc{L) as a 
function of the lattice size L = = are listed in Table [H We plot them 
in Fig. m where the plusses refer to the Polyakov line averaged along the extra 
dimension, Eq. (12. 6p and Eq. (12. 7p . and the crosses refer to the Polyakov line 
evaluated in the first slice along the extra dimension, Eq. (12. 8 p and Eq. (12. 9p . In 
order to define the error of the derived observables Pc{L) and Xc{L) at the critical 
point we used again the bootstrap procedure. We studied both the bootstrap 
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|3„PolyT(X5=lS . 
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Figure 4: Finite size scaling analysis for the transition at 7 = 1, Lj- = 2 based on 
Tabled! We use two definitions of the temporal Polyakov loop, Eq. (I2.6P (plusses) and 
Eq. p.Sp (crosses). The lines are fits to the data using the critical exponents of the 
four-dimensional Ising model. 




1.475 1.48 1.485 1.49 1.495 1.5 1.505 



Figure 5: Finite size scaling analysis for the transition at 7 = 1, Lt = 2: the Binder 
cumulant i?4 of Poly 7^ as a function of /3 for several values oi L = = L^. 

analysis of the reweighted susceptibility and the distribution of the critical value 
fitted on each bootstrap sample. The first analysis is used only to define a fitting 
range for the second one. The critical behavior of a second order phase transition 
is 

Xc{L) = XPoiy.mL)) ~ (3.2) 
|/3,(L)-/3,(L = oo)| ~ (3.3) 

We do not have enough data to predict in a reliable way the critical exponents. 
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Idata li,=P5=1.7, L^=2| 

») 4d string 
») 5d string 

0.6 0.7 0.8 0.9 1 

r/r„ 

Figure 6: The coupling aqq and the slope c on a isotropic lattice with L5 = 2 at /3 = 1.7. 

But what we can clearly state is that all our predictions are totally compatible 
with a phase transition of second order and the scaling behavior dictated by 
the critical exponents of the four-dimensional Ising model, where 7 and u are 
respectively 1 and 1/2, as expected from [H]. The lines in Fig. [Hare fits to our 
data using the critical exponent of the four- dimensional Ising model. For we 
fit the data for L > 4 and the fits work well. For Xc we fit L > 6 but there are 
still some non-scaling features which means that we would need bigger volumes. 
We estimate Pc{L = 00) = 1.48280(1). There is perfect agreement between the 
determination of Pc using both definitions of Poly^-, see Table [TJ We emphasize 
that the two definitions of the Polyakov line must (and indeed they do) lead to the 
same value for the critical coupling in infinite volume and the same value for the 
critical exponents. We notice that the deviations from scaling are larger for the 
definition without the average along the extra dimension and the autocorrelation 
times are bigger but still under control. 

In Fig. [5] we show the Binder cumulant B4 of the temporal Polyakov loop 
Eq. (12.1 op (here we use only the definition of Poly^ in Eq. (12. 6p ). We plot the 
directly simulated (symbols in the legend) and the reweighted data together. The 
Binder cumulants tend to intersect at the value f3c{L = 00) marked by a vertical 
line. Looking at Fig. on the right hand side of the intersection "point" the 
smallest volume [L = 4) corresponds to the lowermost data points and the largest 
volume (L = 16) to the uppermost data points. The value of the intersection is 
for our small lattice sizes L still far from the analytic estimate for the universality 
class of the four-dimensional Ising model 0.27 (horizontal line in Fig. [B], 
indicating again that larger values of L would be necessary for a complete analysis. 
As discussed in [6] the values of B4 computed at (3c{L = 00) approach B^ with 
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corrections which slowly decrease like 1/ ln(L). For our volumes an extrapolation 
of S4 at L = oo is not feasible. Our strongest evidence for the universality class 
of the four-dimensional Ising model comes from the scaling analysis presented in 
Fig.H 

We have computed the static potential V{r) on 32^ x 2 lattices at (3 = 1.7. 
We ran 4 replica for a total of 4 x 10^ measurements. Each replicum used 64 cores 
of the supercomputer Cheops of the University of Cologne and consists of 10^ 
update iterations (each iteration does one heatbath and 16 overrelaxation sweeps) 
for thermalization and 10^ measurements of the Wilson loops separated by 10 
update iterations. We take two levels of HYP smearing for the spatial links of 
the Wilson loops. We measure the potential starting from distance r = 2a and 
the force from r = 2.5a. For the fit of the effective masses in Eq. f l2.14p we use 
the range t = 2, 3, . . . , 6. We obtain the scale ro/a = 5.21(7) with an integrated 
autocorrelation time rint(T'o) = 0.5. The potential can be fitted for r > ro by 
the form predicted from the effective bosonic string theory and we estimate the 
string tension to be oTq = 1.6(3). The left plot in Fig. [H] shows our data for 
the coupling a;qq(l/r) defined in Eq. fl2.17p as a function of r/vQ. For comparison 
we plot the 2-loop (dashed) and 3-loop (solid) perturbative curves for the SU(2) 
Yang-Mills theory in four dimensions, as explained in Appendix |Al The data at 
our smallest distance agree with the 3-loop perturbative curve. The results for 
the slope c(r) defined in Eq. fl2.18p are shown in the right plot of Fig. El There 
is a trend towards the value c(oo) = — 7r/12 predicted by the effective bosonic 
string in four dimension (marked by the dotted line), but the statistical errors are 
already too large at distance r = tq. 

4 Anisotropic couplings 

In order to outline the phase diagram of the system at 7 7^ 1 we start from the 
already known results for 7 = 1 and follow the evolution of the various "criti- 
cal points" as function of /34 and for the different lattice geometries that we 
presented. The critical lines can be cataloged in the following groups: 

• bulk phase transitions. 

This kind of transition seems to be always present if the geometry of the 
lattice is large enough, i.e. in the "infinite volume" and zero temperature 
limit. More precisely for any fixed 7 we find a window of values of /34 and /S^ 
in which we detect a bulk phase transition signaled by an hysteresis effect in 
the quantities Plaq4 and Plaqg. The hysteresis is seen provided the volume 
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/34c 


/35c 




ro/a4 


1.865 


1.34 


10 




1.955 


1.24 


10 




2.33 


0.937 


14 


4.774(44) 


2.5 


0.8697 


20 


8.46(11) 



Table 2: Critical couplings for the bulk phase transition and estimates Lg for the 
minimal lattice size L™™ required to be in large four-dimensional volume and see the 
first order transition. Where we determined it, we list the values of the scale ro/a4. 

is large enough!^: 

Ls > L5mm(7) and min{Ls , Lt} > Lf''{l) ■ (4.1) 

• Center breaking phase transitions in the temporal direction (or analogously 
in any other space direction other than the extra dimension). 

This kind of transition can be obtained by considering a small temporal 
direction Lt < L™™(7) and the other directions larger than their respective 
minimal sizes. The order parameter is the Polyakov loop Poly^. 

• Center breaking phase transitions in the extra direction. 

This kind of transition is obtained with a small extra dimension L5 < 
L5mm(7) while the other dimensions are larger than L™™(7). The order 
parameter is the Polyakov loop winding along the M = 5 direction. 

4.1 Bulk phase transitions 

Fig. [7] shows the behavior of the bulk phase transition in the {^4,(35) parameter 
plane. The points (plusses) have been determined by simulations and their "er- 
rors" are the width of the hysteresis on the largest lattices we have simulated and 
are therefore only indicative. The phase transition is everywhere first order and 
we draw a grey band through the points to guide the eye. We put values for Lg 
(that are valid also for Lt) and L5 which are sufficiently large to see the hysteresis 
signal in the plaquettes and are therefore our estimates for the minimal sizes i^™'" 
and Lbmin. The dotted line simply represents the couplings corresponding to 
7 = 1. 

In [12] it was claimed that at P4 = 3.0 and = 0.779(1) the bulk phase 
transition turns second order, as it was found in a mean-field computation [T5| 

^ We extend the definition of the minimal size Eq. in order to keep track of the different 
lattice spacings and 05. 



15 



4.5 



\ Bulk Phase Transition 

4 \ 



3.5 - 
3 - 

cM? 2.5 - 
2 - 
1.5 - 



0.5 



L^4 



L=8 "-^ L^L=6 

^ U=10 



L3=14 

1 - L =4 , ^ - 

5 + 

L=2, 



0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 

P4 



Figure 7: Bulk phase transitions (plusses, the grey band is to guide the eye) of the 
five-dimensional SU(2) gauge theory using the anisotropic Wilson plaquette action. We 
indicate estimates Lg and L5 for the minimal values and L5min, see Eq. (j4.ip . 

required to see that the transitions in this plot are first order. 



[16]. At /34 = 3.0 lattices of much larger size than the ones simulated so far are 
needed in order to determine the order of the phase transition. Given the large 
computational effort, at present we cannot comment on the situation at (S^ = 3.0. 
In Fig. [H]we present our results at (3^ = 2.5. We show the behavior of Plaq4 as 
a function of for simulations with a hot start or a cold start. Each simulation 
consists of 10"^ measurements separated by 10 update iterations, composed by 1 
heatbath and Ls/2 overrelaxation sweeps each. We keep = and L5 = Ls/2. 
Fig. IH] shows that for Lg = 12 (asterisks for the hot and empty squares for the 
cold run) there is no hysteresis but a smooth cross-over. For = 20 (plusses 
for the hot and crosses for the cold run) the hysteresis appears indicating the 
presence of a bulk first order phase transition. We list some values of the critical 
couplings for the bulk phase transition at 7 < 1 in Table |2j These couplings are 
chosen to be approximately in the middle of the hysteresis region. Table [2] also 
contains the lattice size Lg which is sufficiently large to see the hysteresis and is 
our estimate of L™™ in Eq. (14. ip . If Lt or Lg is chosen smaller than L™™ one 
sees a cross-over which is due to "compactification" of the temporal or spatial 
dimension respectively. 
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Figure 8: Scan in (3^ of the plaquette Plaq4 at /34 = 2.5. Lattices with Lg and Lt equal 
to 20 or larger are needed to see the hysteresis. 

4.2 Center breaking phase transitions 

In Fig. |9]we complete our phase diagram of five-dimensional SU(2) gauge theories 
in the (/34,/95) coupling plane of the anisotropic Wilson plaquette gauge action. 
In addition to the first order bulk phase transitions (empty squares), studied in 
Section H?T| we plot phase transitions due to center breaking (or compactification) 
along either the temporal or the fifth dimension. The order parameter is the 
Polyakov loop in the temporal or fifth dimension respectively. Our new results 
are mainly at 7 < 1 (i.e. /34 > /Ss, the region below the dashed line in Fig. [9]). 
For 7 > 1 (i.e. < (3^, the region above the dashed line in Fig. ^ we include 
the points of [3j. In [21EIES] the transitions due to compactification of the fifth 
dimension have been studied and found to be second order. These transitions are 
represented by empty circles in Fig. [91 For L5 = 2 they extend the 7 = 1 transition 
which we discussed in Section [31 This line ends at 7 < 1 around /34 = 2 when it 
"hits" the bulk phase transition line. The transitions at L5 = 4 and L5 = 6 found 
in |3] branch off the bulk phase transition line at 7 > 1 as /34 is lowered more and 
more. 

In this article we study the transitions due to compactification of the temporal 
dimension. For Lt = 2 (filled circles in Fig. [9]) they are again continuation of the 
7 = 1 transition. We performed a finite size scaling analysis for the transition 
point at /Ss = 0.5. In Table [21 we list the critical values /34c of the coupling f]^ 
at which the susceptibility of the temporal Polyakov line has its maximum Xc for 
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Figure 9: The phase diagram of the five-dimensional SU(2) gauge theory with the 
anisotropic Wilson plaquette action. Explanations are in the text. 

several values of the lattice size Lg = 8 . . . 24 keeping L5 = Ls/2. We use both 
definitions of the Polyakov line averaged along the extra dimension Eq. (12. 6 p or 
taken on a fixed slice X5 = 1 Eq. ( 12. 8p . In Fig. [10] we plot the data as function of 
l/Lg together with fits using the critical exponents of the four-dimensional Ising 
model, cf. Eq. fl3.2p and Eq. f l3.3p with /3c replaced by and L hj Lg. The fits 
work quite well and using the largest three volumes we estimate the critical value 
(34c{Ls = 00) = 1.82115(8). Both definitions of the Polyakov line give a perfectly 
compatible result. 

In Fig. [TT] we show the results for the Binder cumulant i?4 of the temporal 
Polyakov loop Eq. (12.101) (here we use only the definition of Poly^ in Eq. (12. 6p ). 
We apply the reweighting technique as discussed in Section [3] to get a denser data 
set. There is clear trend in the data to intersect at a common value. The vertical 
line marks /34c and the horizontal line the expected value for the universality 
class of the four-dimensional Ising model [B]. In order to precisely determine the 
critical behavior larger volumes are needed. Please note that the lattices which 
we simulated already require significant computational resources. 

For Lt = 4: (filled squares in Fig. |9]) we faced a new situation. Staying at 
/^s = 0.5 (like for L-r = 2), in the broken phase close to the phase transition the 
histogram of the temporal Polyakov loop Eq. (I2.6p shows several peaks. This is 
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Volume 


/34c (Poly^^) 


/34c (PolyT(l)) 


Xc (Poly^) 


Xc (Polyr(l)) 


2 X 8^ X 4 
2 X 12^ X 6 
2 X 16^ X 8 
2 X 20^ X 10 
2 X 24^ X 12 


1.83535(14) 
1.82713(19) 
1.82455(20) 
1.82325(7) 
1.82263(3) 


1.83751(25) 
1.82859(20) 
1.82588(19) 
1.82393(8) 
1.82309(4) 


46.38(18) 
112.6(1.4) 

216(5) 
351(5) 
525(6) 


16.69(7) 

29.3(3) 

44.4(9) 

57.4(6) 

70.9(5) 


2 X oo^ 


1.82115(8) 


1.82095(19) 







Table 3: For the center breaking phase transition at (3^ = 0.5, = 2, we Hst the 
critical values of the coupling /34c and the susceptibility of the temporal Polyakov loop 
Xc as function of the volume. We use two definitions of the Polyakov line, Eq. (j2.6p and 
Eq. dm). 



Volume 


/34c (Polyr(l)) 


Xc (Polyr(l)) 


4 X 123 X 6 


2.3024(8) 


23.44(16) 


4 X 16^ X 8 


2.2993(5) 


41.1(3) 


4 X 20^ X 10 


2.2995(6) 


62.2(5) 


4 X 24^ X 12 


2.2985(10) 


88(3) 


4 X 32^ X 16 


2.2981(4) 


156(2) 


4 X oo^ 


2.2974(4) 





Table 4: For the center breaking phase transition at (3^ = 0.5, Lt = 4, we list the 
critical values of the coupling /34c and the susceptibility of the temporal Polyakov loop 
Xc as function of the volume. We only use the definition of the Polyakov line Eq. (j2.8p . 

shown by the plot on the left of Fig. [12] for /34 = 2.32. If we do not average the 
Polyakov loop along the fifth dimension and take its definition Eq. (12. 8 p for the 
slice X5 = 1, we get the histogram shown by the plot on the right of Fig. [121 There 
is only a single peak. The interpretation of these plots is that we are in a phase 
of the theory where the Polyakov lines can fluctuate almost independently in the 
hyperplanes orthogonal to the fifth dimension. In each hyperplane the distribution 
of their values has the characteristic double-peak shape which appears like the plot 
on the right of Fig. [T2] when taking the absolute value. When averaging over the 
hyperplanes they produce a multiple-peak structure, which is therefore an artefact 
and does not mean the presence of multiple (more than two) vacua. This is a very 
interesting effect signaling that the interactions between the hyperplanes is weak. 
The decoupling of the hyperplanes for 7 < 1 was in fact also found by the mean- 
field computation of [TB]. While the geometrical setup of our computation is 
different, the qualitative feature we see are the same. 



19 



he Polyr ' 



Xpolyj ' 
XpoIi.t(Xs.1) 



300 
200 



0.05 0.06 0.07 0.08 0.09 0.1 0.11 0.12 0.13 
1/L 



0.04 0.05 0.06 0.07 0.08 0.09 0.1 0.11 0.12 0.13 

1/U 



Figure 10: Finite size scaling analysis for the transition at (3^ = 0.5, Lt = 2 based on 
Table O We use two definitions of the temporal Polyakov loop, Eq. (I2.6P (plusses) and 
Eq. p.Sp (crosses). The lines are fits to the data using the critical exponents of the 
four-dimensional Ising model. 




Figure 11: Finite size scaling analysis for the transition at /^s = 0.5, Lt = 2: the Binder 
cumulant B4 of Poly^ as a function of /34 for several volumes. 

Taking the temporal Polyakov line defined only on the slice X5 = 1 we can 
do a finite size scaling analysis as we previously did for Lt = 2. It is now clear 
why we defined two Polyakov loop operators. Indeed in the present case only for 
the one not averaged over the extra dimension a finite size scaling analysis can 
be made. Where it was possible to perform both analyses, we showed that their 
results coincided. The critical values /S^c and the maximum of the susceptibility 
Xc are listed in Table H] and shown in Fig. [131 The fits to the critical behavior 
Eq. (13. 2p and Eq. (13. 3p assuming the critical exponents of the four-dimensional 
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Figure 12: Histograms of the absolute value of the temporal Polyakov loop at /34 = 2.32, 
/Ss = 0.5 and Lt = 4: taking the definition Eq. ()2.6p (average over the extra dimension, 
left plot) or Eq. (j2.8p (at a fixed slice X5 = 1, right plot). 
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Figure 13: Finite size scaling analysis for the transition at (3^ = 0.5, Lt = 4 based on 
Table HI The temporal Polyakov line is taken in the slice X5 = 1 . The lines are fits to 
the data using the critical exponents of the four-dimensional Ising model. 

Ising model work well and give (5ac{.Ls = 00) = 2.2974(4). 

In summary the phase transitions due to compactification of the temporal 
dimension at 7 < 1 for L-r = 2 and Lt = 4 are second order and compatible with 
the critical exponents of the four- dimensional Ising model thus confirming [H] . 

4.3 The static potential 

The simulation results presented in the previous section show that at 7 < 1 
dimensional reduction from five to four dimensions can occur by compactifying one 
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Figure 14: The coupling Oqq Eq. (j2.17p at the bulk phase transition in the confined 
vacuum. The simulations have been performed at /34 = 2.33, 2.5, cf Table [2j 



of the directions orthogonal to the fifth dimension. We discovered that the relevant 
degrees of freedom (the broken Polyakov loops) fiuctuate almost independently 
M^hen defined in the hyperplanes orthogonal to the extra dimension. The theory 
in the broken phase reduces to four dimensions and not three, thereby indicating 
that the hyperplanes do not exactly decouple but some interaction is left. Here 
we go back to simulations of the theory in "infinite volume". By this we mean 
that all five directions are larger than their minimal size, see Eq. f l4.ip . If the 
hyperplanes decouple, the theory reduces to four dimensions, as is the case in the 
mean-field computation of [T6] . 

We computed the static potential V{r) along spatial directions orthogonal to 
the extra dimension (and averaged over the extra dimension). In the deconfined 
phase (where the Polyakov loops are broken) it is a five-dimensional Coulomb 
potential [13] • Here we choose parameters which correspond to bulk phase tran- 
sition points and we put ourselves approximately in the middle of the hysteresis 
curve. We started the simulations with a hot start in order to stay in the vac- 
uum of the confined phase and we checked that the Polyakov loop expectation 
values are indeed zero in all directions. We choose the points in parameter space 
(/34 = 2.33, /35 = 0.937) and {P^ = 2.5, /Jg = 0.8697), cf Fig. El and the lattice 
size is 32^ x 16. We ran 4 replica at (S^ = 2.33 and 11 replica at = 2.5 for 
a total of 38280 and 103410 measurements respectively. Each replicum used 512 
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Figure 15: The slope c Eq. (|2.18p at the bulk phase transition in the confined vacuum, 
from the same simulations as in Fig. [TH 

cores of the supercomputer Cheops of the University of Cologne and consists of 
10^ update iterations (each iteration does one heatbath and 16 overrelaxation 
sweeps) for thermalization and approximately 10^ measurements of the Wilson 
loops separated by 10 update iterations. We take four levels of HYP smearing 
for the spatial links of the Wilson loops. We measure the potential starting from 
distance r = 2a4 and the force from r = 2.5a4. For the fit of the effective masses 
in Eq. fl2.14p we use the range t = 2, 3, . . . , 6. The values of the scale VQ/a^ are 
given in Table O At = 2.5 the lattice spacing is almost half the one at 
/3i = 2.33. 

In Fig. [H] we show the coupling aqq Eq. f l2.17p as a function of r/ro, square 
symbols for = 2.33 and circles for P4 = 2.5. The string tension is estimated 
to be ar^ = 1.3(2) and otq = 1.4(2) at /34 = 2.33 and P^ = 2.5 respectively. 
We plot the 2-loop (dashed) and 3-loop (solid) perturbative curves for the SU(2) 
Yang-Mills theory in four dimensions, as explained in Appendix [Aj The points at 
the smallest distances are compatible with the 3-loop running, especially on the 
finer lattice at = 2.5. 

Fig. [T^ presents the analysis of the slope c(r) Eq. fl2.18p . At P^ = 2.33 the 
values of c seem to be more compatible with the five- dimensional value of the 
Liischer coefficient (marked by a horizontal dashed-dotted line). As we lower the 
lattice spacing 04, at (3^ = 2.5 the data have a clear trend towards the four- 
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dimensional value of the Liischer coefficient (marked by a horizontal dotted line). 
We loose the statistical signal above distance r = tq. But as was noticed in [35] 
the onset of the string behavior in the static potential starts already at distance 
tq. There is no scaling between the two simulations and we are most probably 
not on a line of constant physics. We interpret our results as an indication that 
increasing /34 along the bulk phase transition line makes the theory more four- 
dimensional than five-dimensional, hinting at a dynamical localization mechanism 
for the gauge field. 

It is not clear whether there exists a continuum limit when /34 increases while 
staying in the confined phase at 7 < 1. But even if it does not exist, there 
could be a window of values of the lattice spacing for which the cut-off effects 
are small and the five- dimensional gauge theory can be described by an effective 
four-dimensional theory for energies much smaller than the cut-off. 

5 Conclusions 

Lattice gauge theory is the tool to explore five- dimensional gauge theories away 
from their trivial limit. So far we did not find in Monte Carlo simulations of gauge 
group SU(2) using the anisotropic Wilson gauge action and periodic boundary 
conditions a second order bulk phase transition where a continuum limit could 
be taken. We located a line of ffist order bulk phase transitions separating the 
confined from the deconfined phase. In this article we studied in particular the 
phase diagram for anisotropy 7 < 1 (where the lattice spacing 05 is larger than 
04). As we move along the bulk phase transition line, while decreasing 04 and 
staying in the large volume confined phase, we find indications of dimensional 
reduction. Dimensional reduction at 7 < 1 was previously found in a mean- 
field calculation [16] and relied on a decoupling of the hyperplanes orthogonal to 
the extra dimension. This effect is seen in our analyses of second order phase 
transitions related to breaking of the center. At these transitions our finite size 
scaling studies are perfectly compatible with the critical exponent of the four- 
dimensional Ising model. 

As a next step we would like to measure the mass spectrum in the scalar 
and vector channels in order to understand what exactly is the (almost) four- 
dimensional theory that we get in the hyperplanes orthogonal to the extra dimen- 
sion at 7 < 1. These studies prepare the ground for future simulations of the 
theory with orbifold boundary conditions [2III111I15]- 
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A Perturbation theory in the qq-scheme 

We summarize the steps needed to compute the 3-loop running of 

Att 

= ^r^FW, /i = l/r, (A.l) 

where Cf = (A^^ — 1)/(2A^), in the SU(N) Yang-Mills theory in four dimensions. 
The beta function is defined by the renormalization group equation 

/^^^qq(/^) = P<n{9<iM) (A.2) 
and has the perturbative expansion 

M9^^) = -9'(bl^''^+b'^'^9l + b^2''^9'--)- (A.3) 



^qq^^yqq/ ~ yqq "i" "i yqq "r "2 yqq 

The coefficients (Ca = A^) 



&i'^'^ = &o = (^11^a/3, (A.4) 
6i^^) = 6i = ^34Cl/3, (A.5) 



are universal. The 3-loop coefficient 

r{qq) _ Cl f206 447r2 llvr^ 242 



TT 



^ (47r)6 \ 3 3 12 9^^7V3 J k J 

is determined combining results from [IHIHT]. In order to integrate Eq. f lA.2p we 
evaluate numerically 



(A.7) 

where for /3qq we insert the truncated perturbative expansion. 
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For the evaluation of Eq. ( 1A.7I) we need to know the Lambda-parameter in 
the qq-scheme. For gauge group SU(2) it can be estimated from the data of the 
Schrodinger Functional (SF) coupling in Table 5 of [H]. We integrate Eq. (IA.7P 
in the SF scheme using the coefficient determined in |l9l[50]. We take the 
smallest coupling from |38] and use the determination of the tq scale in [M] . 
After converting to the qq-scheme jlH] we estimate 

Aqqro = 0.80(11) (A.8) 

from the 3- loop running (which agrees with the 2-loop result). 
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